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Abstract 

A version of foliated spacetime is constructed in which the spatial geometry 
is described as a time-dependent noncommutative geometry. The ADM version 
of the gravitational action is expressed in terms of these variables. It is shown 
that the vector constraint is obtained without the need for an extraneous shift 
vector in the action. 

Introduction 

The problem of divergences in quantum field theory, and especially in quantum grav- 
ity, strongly suggests a need to describe the geometry of space-time as something 
different or more general than a classical manifold of points. One interesting general- 
ization of classical geometry is noncommutative differential geometry. It has attracted 
some attention in physics in recent years, mostly in the form of Connes-Lott and re- 
lated models of particle physics . It has also been used as a regularization technique 
for Euclidean quantum field theory 0. Attention to applying noncommutative ge- 
ometry to gravitation has been limited thus far; an expression for the Euclidean 
gravitational action was proposed by Connes and calculated explicitly by Kalau and 
Walze in and by Kastler in ||. 

In noncommutative differential geometry (see []3], ^, |5|, |J) the geometric informa- 
tion is set in the form of a real spectral triple (A,7i, D). This consists of a *-algebra 
(of "functions") A acting on a Hilbert space 7i and an unbounded, self-adjoint Dirac 
operator D — also acting on 7i; the real structure is an antiunitary operator J, and 
for even dimension there is a Z 2 -grading operator on 7i. 



An ordinary Riemannian spin geometry can be described by a spectral triple. 
In that case A is the commutative algebra C°°{M.) of smooth C- valued functions 
with the involution (the *) being complex conjugation; the Hilbert space is the TC = 
L 2 (A4, SM) of spinor-functions with the canonical inner product (tp | tp) — j M <pipe; 
the Dirac operator is the usual D = ry^'V,-; the real structure J is the charge con- 
jugation operator; and the grading is the chiral grading — the generalization of the 
familiar 75. For the details of the axioms a spectral triple must satisfy in general, see 
§■ 

If a real spectral triple is commutative (i. e., the algebra is commutative) and 
satisfies a slightly stronger axiom (that J a J" 1 = a*), then up to equivalence it is the 
spectral triple of an ordinary manifold, although with some unwanted additional free- 
dom to the Dirac operator (see section |). In particular, the point set can be recovered 
as the spectrum of the algebra (as in Gelfand's theorem). In the noncommutative 
case, the spectral triple becomes a generalization of a Riemannian spin manifold. 

We would like to use noncommutative geometry to generalise the geometry of 
space-time, but there is one rather obvious obstacle to this. Space-time at macroscopic 
scales resembles a Lorentzian rather than Riemannian manifold, and noncommutative 
geometry cannot straightforwardly describe a Lorentzian space-time. There is no 
natural positive definite inner product for spinor-functions on space-time, and if an 
inner product is chosen anyway the Dirac operator cannot be self-adjoint since it has 
complex eigenvalues. 

There are many techniques of analysis which do not really work for space-times - 
enough that there is a standard way of makeshifting around the problem. The space- 
time is foliated into space-like hypersurfaces, and the leaf space is parameterised by 
a time variable. Since the leaves are spacelike, the geometric situation is reduced to 
Riemannian geometry. 

The standard Hamiltonian approach to general relativity is the ADM (Arnowitt- 
Deser-Misner) formalism (see f|TTJ ) . The intrinsic geometry of the leaves is simply 
described by a Riemannian metric. Two types of information describe how the leaves 
fit together. The lapse function gives the infinitesimal proper time separating the 
leaves at each point. The shift vector is needed because, rather than treating the 
leaves as separate topological spaces, the Hamiltonian approach treats them as a 
single manifold with time-dependent geometry. This means that the separate leaves 
must be identified pointwise, and the identified points trace out curves in space- 
time. The shift vector is the velocity of these curves, measuring their deviation from 
orthogonality to the foliation; it parameterises the arbitrary choice made in identifying 
points. By this construction, a time-dependent Riemannian metric, scalar field, and 
vector field on an n-dimensional manifold determine an n + 1-dimensional space-time 
geometry. 

The construction described in this paper is an algebraic Hamiltonian formalism 
- a combination of noncommutative geometry and the standard Hamiltonian for- 
malism. The time-dependent Riemannian geometry can obviously be described by a 
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time-dependent, commutative spectral triple (At,Ti.t,T> t )- The spacing of the leaves 
continues to be described by the lapse which is now thought of as taking values in At 
at each time. In the ADM formalism, the variables were time-dependent functions 
and fields, and it was necessary to fix the manifold on which these live. In this case 
the variables are operators and it is necessary to fix the Hilbert space on which these 
act. This means identifying all the TCtS to a single TC. One way to do this would 
be to extend the pointwise identification to identify the entire spinor bundles of the 
different leaves. This would be essentially a Schrodinger representation, and would 
make the algebras At all identical. 

A much more covariant approach is to take the Heisenberg representation. The 
spinor-functions at different times can be naturally identified by the 0-mass space- 
time Dirac equation. This identification involves no arbitrary choice and thus requires 
nothing like a shift vector. 

In the ADM formalism the action is given by a Lagrangian which is equivalent to 
the Einstein-Hilbert Lagrangian modulo a classically irrelevant overall time derivative. 
The ADM Lagrangian can be split into kinetic and potential terms. The potential 
term is proportional to the integral of the spatial scalar curvature weighted by the 
lapse; an algebraic expression for this is easily adapted from Connes' expression for 
the Euclidean gravitational action. The kinetic term can be expressed in terms of the 
lapse and the extrinsic curvature tensor. The extrinsic curvature is simply related 
(via the lapse and shift) to the time derivative of the metric tensor, so in some 
sense it comes from the time derivative of the intrinsic geometry. In the algebraic 
formulation the intrinsic geometry is described by the Dirac operator; in order to 
extract the information of the extrinsic curvature which is needed to calculate the 
Lagrangian, we should look to the time derivative of the Dirac operator. 

In section [I] most of the basic structures, notations, and formulas are introduced. 
In |2| it is described how to calculate time derivatives of operators, and this is applied 
to the Dirac operator in section |3|. In |] the construction of the real structure and 
gradings are discussed as well as the issue of what type of spinors are used. The 
noncommutative analogue of integration is described in section |5|. In section |6| the 
action and dynamics of a scalar field are described in this geometrical setup; this is a 
simpler case of some of the issues that come up with the gravitational action which 
is calculated in section [7]. The nature of diffeomorphism invariance is described in 
section || and it is seen how the vector constraint comes about in this model. The 
remaining issues necessary to generalize this to truely noncommutative geometry are 
discussed in section || The failure of time diffeomorphism invariance in this model is 
discussed in section [H]. The qualities and possible use of this model are discussed in 



section 11. The appendix describes a simple calculation of the overall normalization 



of the Wodzicki residue and the normalization of the Euclidean gravitational action. 
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1 Basic Definitions 



Let M = K x S be a space-time such that each of the hypersurfaces Et := {£} x £ is 
spacelike (in fact, Cauchy) and compact. These hypersurfaces will be regarded as the 
leaves of a foliation. Assume that M. has a spin structure. The signature convention 
is time-like. Greek indices are space-time indices; Latin indices are spatial indices 
( "space" = E 4 ). When needed, an " M " will distinguish a space-time object from a 
spatial object. V and | will denote the space-time covariant derivative; D and ; will 
denote the space covariant derivative. A slash denotes a tensor contracted with the 
Dirac vector 7 M , as in gt = 7 M a M . 

For any two spinor-functions if, if) G C°°(S t , SM) there is a positive definite inner 
product given by the flux of the current (p *7"0; this generates a Hilbert space 7i t — 
L 2 (H t , SAi) for each t. For solutions of the 0-mass Dirac equation = this 
current is conserved, and so the inner product is independent of the specific S t at 
which the flux is taken. This means that solutions of the Dirac equation span another 
Hilbert space 7i. For every t, there is an obvious map Ti — > 7i t given by restricting 
to E t . This is actually an isomorphism; an arbitrary spinor-function at £ f is valid 
initial data for the Dirac equation, so propagation gives the inverse map 7i t ~^ 
These isomorphisms naturally identify each 7i t with 7i. 

Each of the algebras At '■— has an obvious action on 7i t and hence on 

7i. Ar will be regarded as a bundle of algebras over R. Since the E t 's are all 
diffeomorphic, the algebras At are all (non- naturally) isomorphic. 

The first calculations here are most conveniently done on Dirac spinors since some 
of the operators involved reverse chirality, but I will soon restrict to chiral spinors 
(section ^). For odd spatial dimension (hence even space-time dimension) there is a 
chiral grading r of 7i; this is the generalization of the usual 75. For all dimensions, 
there is an antiunitary charge conjugation operator C acting on 7i. Both T and C are 
time-independent, since they preserve solutions of the 0-mass Dirac equation. 

The inner product can be written {(p\ifj) = J St <p *^ifj = J St ftftipc, where n is the 
unit normal to E t and e = n j M e is its volume element. This means, in particular, 
that for any Dirac matrix valued function A, the adjoint is A* = jlAyl; so, ji* = fl 
and 7; = -7j. 

Let a J be the Dirac vector of E t . The spatial Dirac operator is T> := ia^ Dj = 
ia^(dj + ...). It should be the Hamiltonian associated to some sort of normal 
derivative; that is, the Dirac equation can be written in a Hamiltonian form as 
i(d n + . . = Vip. On the other hand = ^tp = (fid n + j j dj + . . .)ip, so 
id n ip = (—iflj^dj + . . . )ip; therefore, a J = — fi^ = ^ yi = (« J )*. This satisfies 
the required identity [a^ , a k ]+ = 2g^ k with the positive definite metric of E t . Analo- 
gously to the slash notation, a tensor contracted with a will be denoted with a hat. 
Certain spatial tensors will also be used as space-time tensors, the identification is 
made in the "natural" index positioning of the tensor. 

The V derivative of the normal vector is rz/f = K^. + n„a M . K^. is the extrinsic 
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curvature of the leaves (in natural position), and a M := n^ n is the acceleration of the 
unit normal (also in natural position). Both of these are really spatial tensors - 
they vanish when contracted with n — and are written here in natural position. In 
covariant positon is symmetric. Using a differential form notation, the lapse N 
is the scalar defined by the equation n = Ndt; since a is spatial, we have 

a = n j (n A a) = n jdn = n j d(Ndt) = n j (dN A dt) 
= — 7V _1 n j(nA dN), 

i. e., the projection of —N~ x dN orthogonal to n. Since a is naturally contravariant 
rather than covariant, we have a? = N~ 1 N' : ' . 

Parallel transport by D respects the metric and keeps tangent vectors tangent to 
the foliation, therefore D regards the unit normal vector as constant; so = [Dj, ^]_. 
In general, A i— > \$[yi, A] + is a projection onto operators commuting with yi (taking 
advantage of fl 2 = 1), and Dj (acting on spinors) is exactly the image of V,- under 
this projection: 

Dj = \i[Vj, i\ + = V, + \i[Vj, i\- = v, + \ii^ 3 

= v 3 + \k r (i) 

Kj := OL^K k j is naturally covariant and self-adjoint. Consistency with the inner 
product means that Dj must be skew-adjoint; in fact —\Kj and Dj are precisely 
the self-adjoint and skew-adjoint parts of V,-. 
The spatial Dirac operator is 

V = ia j Dj = ia ] (Vj + \K^j = -ijt*fVj + \K 

= -liy + lV n + \K (2) 

(where K = tr K). It is self-adjoint, anticommutes with ^, and commutes with T. 

2 Time Derivatives 

In order to find an expression for the gravitational action, we need to be able to 
extract the extrinsic curvature from the time derivative of V. To do this we need an 
expression for the time derivative of T> in terms of more familiar objects. For this 
we need a way of calculating time derivatives of operators. Toward this end, first 
consider spinor-functions on M.. 

A space-time spinor-function ip E C°°(M, SM) restricts to a spatial spinor- 
function T, t — > SM. for every value of t E R; using the identification of H t with 
7i, this defines a function R — > 7i (i. e., t h- > i(}(t)) which is completely equivalent to 
ip. One might ask what operator on C°°(M., SA4) corresponds to ^ on C°°(R, TL). If 
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if G C°°(Ai, SM) is a solution of the Dirac equation (that is, ip G H), then the inner 
product with ip taken at S t is the inner product with ip(t), i. e., {f\ip)t = (fli^it))- 
We need to compute 4r{ip\ip) t = (ip\-^ip(t)). To do this, use Stokes theorem on the 
infinitesimal volume between E 4 and £t+dt, and the formula M e = Ndt A e, 

= / tpNflfflipe 

JT, t 

= (cp\Nf[yip) t . 

So | = Ntf. 

An operator A that acts on space-time spinors can be applied to a solution ip of 
the Dirac equation to give a space-time spinor Aip. With the above identification Aip 
can be thought of as a time-dependent element of 7i. This means that A determines a 
time-dependent operator on 7i. This construction gives a projection from space-time 
operators to time-dependent operators, and it has a preferred left inverse. Given a 
time-dependent operator on H, there is a unique operator in its preimage that only 
involves spatial derivatives. If time-dependent operators are identified with space-time 
operators of this type, then time derivatives are easily calculated; the time derivative 
of an operator is A = [~, A]-. 

In particular, suppose / is a scalar function on M.. Then 

/ = [AW, /]- = JW/|„ = Nf ]n - Naif.,. (3) 

Clearly, this is quite different from choosing a specific time direction and taking the 
partial derivative of /. 



3 The Time Derivative of T> 

Now we can calculate 

v= [Niy,v\_ = -[v,N]_iy + iv[^y,i?]_ 

using the fact [D, iV]_ = iNa, and the formula (0) for T>. Now, = $ and the 
definition a = [V n , n]_ give 

V = -iN0 + iAW + iNfi[f , V n + 
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We can evaluate part of this using the identity that [V M , V^]_ = \ M R Ufl : = \j a 'j l3M Rai3^ 
This gives 

[y,v„]_ = [y,n^]_v M + 7 ^[v„v M ]- 

= [y,^]_v M + i^x ( 4 ) 

since Y M ft^ = rri^Rap^ = ( M 9 va - M g^+^ + yl^-fl) ^R a ^ = 2 M R^ 
(because of the cyclic identity). Now use V • n — K and the definition of M R to get 

which can be inserted back into (|]) to give 

[y,v„]_ = |[v,jy,^]_] + -i[y,K]_, 

and in turn gives 

P=iJV^[V MJ [y,n"]_] + . 

We need to reexpress this in terms of D rather than V; the only difficulty is the 
action of D on the vector index. Note that [y , = fi{a^ + K^) is orthogonal to n 
(hence spatial). If v is a (scalar valued) spatial vector, then 

(Ndt A e)v^ = M e = d(v _j M e) = d{vj [Ndt A e]) 

= dt A d(Nv j e) = dt A e (Nv j ).j 

and so [V M , = i^, = N~ 1 (Nv j )-j. Since by (|) = ^'V,- = ^'(-Dj - this 

gives y^ = N~ 1 (Dj — ^Kj)Nv : ' which continues to be valid when v is matrix-valued. 
This gives 

V = \i\D 3 - \K h Nfi(a j - K j )] + 

= i[D jt N(a j - k>)} + + \N[K V a j - 

since fi commutes with Dj and anticommutes with Kj. The last term vanishes to 
give 

V = \[D v N j - NK j } + = i[Dj, [D\ JV]_ - NK j ]+ 

= i[I?,N]--i[D j ,Nk'] + 

but D 2 = —T) 2 + \R and R — being a scalar — commutes with N. This gives the 
final result: 

V = -\ [V 2 , iV]_ - | [Dj , NK j }+ =: V + + V_. (5) 
These two terms, respectively, commute and anticommute with ft. 
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4 Dirac and Chiral Spinors 



The charge conjugation operator C commutes with the covariant derivative operators, 
and anticommutes with the imaginary unit i. Depending on the dimension, C some- 
times commutes and sometimes anticommutes with 7 M (see fllU|), but this sign cancels 
itself so that C always commutes with ot> . Together this means that C anticommutes 
with D. To construct a real spectral triple describing the spatial geometry, we need 
an antiunitary operator J which should (among other conditions) anticommute with 
D for spatial dimension n = 1 mod 4 and commute with it for other dimensions. 
Since yi anticommutes with T>, this condition is satisfied if we set J = C for n = 1 
mod 4 and J — flC otherwise. This choice also satisfies the other conditions for J, 
and leads to JT = TJ. 

If the space dimension is even, then the Dirac spinors on M. have the same di- 
mension as Dirac spinors on E. In fact the Dirac spinor bundle SM. restricted to T, t 
is equivalent to SH t . 

On the other hand, if the space dimension is odd, SAi restricted to E t is equivalent 
to two copies of SH t . However, chiral spinors on Ai restrict to a single copy of SUf 
In order to work with the Hilbert space of Dirac spinor-functions on E, we must 
restrict H to be the ±1 eigenspace of T (left or right handed spinors). This can be 
done, since all operators of interest from now on are T-even (commute with T). 

If the space dimension is even then there is no T, but a grading operator is needed 
for an even dimensional spectral triple. If we temporarily work in an orthonormal 
frame such that the direction is the normal direction, then the grading is (modulo a 
phase factor) a 1 . . . a n ~ (7°) n 7 1 • • • 7 n - For even dimension n, the product 7 . . . 7™ 
is just a phase; so the grading is ~ (7°) n ~ 1 = 7 = fl, which is precise (up to a sign) 
since it is Hermitian. 

To summarise, we have now expressed the geometry of A4 as a time-dependent, 
real, commutative spectral triple (At,TC, T>t), and a lapse function N(t) G At- If the 
dimension n is even, then the (time-dependent) grading is fi and the real structure 
is J = fiC given by a time-independent antiunitary C. For odd dimension, if n = 1 
mod 4 then J is time-independent, but if n = 3 mod 4 then J does depend on time. 

5 The Wodzicki Residue and Dixmier Trace 

For some vector bundle V — > E over an n-dimensonal manifold, given a choice of 
coordinate system and trivialization, an arbitrary psuedodifferential operator A can 
be written ClS db Fourier integral operator 



in terms of its (total) symbol cr(A) : T*E — > End(l^). Define (Jk{A) as the compo- 
nent of order k in p. The highest order component is called the principal symbol, 
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<7 pT (A) = aord(A)(^4)- Expressing A in terms of its symbol is equivalent to writing 
A as a formal power series in pj := idj with p's always ordered to the right of x's. 
Although the total symbol is strongly dependent upon the coordinates used to define 
it, the principal symbol is actually coordinate invariant; it also has the simple prod- 
uct property a pr (AB) = cr pr (A)a pT (B). This is an example of the pseudodifferential 
calculus; calculations on operators can be performed in terms of their symbols. 

Given an arbitrary Riemannian metric on £, define (<j)(x) as the normalized trace 
of a averaged over the unit sphere in T*T. Because of the invariance of the principal 
symbol, (cr pr (yl)) is coordinate invariant; (<7- n (A))e is also coordinate invariant and 
is independent of the metric. 

The normalized Wodzicki residue is defined as 

jA = J(a- n (A))e. 

It is a trace (i. e., satisfies the cyclic identity) and is proportional to the Dixmier 
trace (see and the appendix), which (for £ compact) can be defined using a 
fiducial Dirac operator D, or any operator such that cr pr |D| = ||p|| (where \D\ : = 
(D*D)^ = (D 2 )^ is the operator absolute value). Define the zeta function Cd.a(s) as 
the analytic continuation of Tr (|Z)|~ s v4). The (extended) Dixmier trace is Tt^(A) = 

Res s=0 ( D n A (s) = ^Res s=0 Cz?,A(s) ( see 0)- 

D can still be used to define the Dixmier trace when D is singular. In general, 
if A is multiplied by a finite codimension projection, then Qd,a will only change by 
a holomorphic function, and Tr w A will be unaffected. It is thus possible to project 
out any finite dimensional eigenspace of D when defining Cd,a, so in particular the 
O-eigenspace of D can be projected out. This should work provided that has finite 
degeneracy in D, as it always does for a compact manifiold. 

The "limiting process" u depends upon the fiducial D, but if A belongs to the ideal 
Cp-'°°)(7i) of operators of differential order —n, then the Dixmier trace is independent 
of the fiducial D and for diagonalizable A is (see 0) 

1 - 

k=l 

where the A's are the eigenvalues arranged in order of descending |A|. The Dixmier 
trace vanishes on the ideal C^' 00 ^ C Z^ 1 ' 00 ) of operators of lower order than — n. One 
of the axioms of noncommutative geometry requires that |X?|~™ e C^ l,co \ 

Now suppose that E is given a Riemannian spin structure, and that V = ST, the 
spinor bundle. The canonical metric and Dirac operator can be used as the fiducial 
ones; this means that \D\ can be used to shift degrees. If we need the integral of 
(a pT (A))e, but A is not of order —n, it is merely necessary to multiply A by an 
appropriate power of \D\. If / G C°°(S) then a(f; x,p) = f(x), so 

jf\D\- n = J fe. 
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Thus a volume integral may be reexpressed in a purely algebraic way, and the Dixmier 
trace may be regarded as algebraic integration. 

In order to evaluate Wodzicki residues, we need to evaluate expressions of the form 
(pi . . . pi). If there is an odd number of p's then this will simply vanish. If there is an 
even number of p's then symmetry implies (piPj . . .pip m ) Qtij ■ ■ ■ 9im)- The propor- 
tionality can be determined by recursion from the fact that p 2 is equivalent to 1 inside 
an (. . . ). In particular (p a p b ) = \g ah and (p a PbPcPd) = (9ab9cd + 9ac9bd + 9ad9bc)- 

The power |D| 2_ri is particularly interesting. From the setup of the calculation in 
terms of symbols, it is readily apparent that cr_„ (\D\ 2 ~ n ) is second order in derivatives 
of the metric (i. e., linear in second derivatives, quadratic in first derivatives) and first 
order in derivatives of the gauge potential (when there is one); it also vanishes for 
flat space with Cartesian coordinates, trivial spin connection, and gauge potential. 
(cr_ n (|-D| 2-n )) is coordinate invariant and gauge invariant, and therefore can only be 
proportional to the scalar curvature R. The precise relation is (see the appendix and 



It is possible to describe matter fields in this geometric context. The simplest case is, 
of course, a scalar field. If is a scalar field on A4, it can be thought of as a function 
4>(t) G A t . Its time derivative is = [f v 4>]- = N^cj)^ = N(j)\ n - Na j (p.j. On the 
other hand [D, 0]_ = ia^(p-j. Combining these gives the Lagrangian for a real scalar 



In the commutative case, there is some freedom to reorder factors in expressions 
of this kind. When an individual term inside a Wodzicki residue is in the ideal £ < - 1 '°°- ) , 
then reordering scalar factors will only change the expression by operators of smaller 
order; these are in £q 1 '°°' ) and vanish when the Wodzicki residue is taken. Functions are 
always of order 0, and in the noncommutative case the commutator of two functions 
is also a function and is thus of the same order; the freedom to reorder is significantly 
reduced. This means that when a Lagrangian is generalized from the commutative 
case to the noncommutative case, there is a factor ordering ambiguity. The factor 
ordering here is chosen as the simplest which makes this Lagrangian manifestly real. 

The canonical momentum conjugate to is ir^ := but what does this mean? 

It is a linear function of the possible variations of <fi for fixed 0. Looking at the above 




(7) 



6 A Scalar Field 



field 
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formula for <j) it is clear that it can only vary by an element of At, so G A* t . It is 
in fact 

Specifying the coefficient of u in this integrand is not equivalent to giving 7r^; the 
coefficient can, for instance, be changed by any operator in Cq , because this will 
not affect the Dixmier trace. This means that ir^ cannot easily be represented by a 
corresponding operator, and there is no simple way to write <fi in terms of and n^. 
This makes a closed form expression for the Hamiltonian essentially impossible. This 
same problem occurs with the gravitational action. 



7 The Gravitational Action 

The principal symbols of the terms of T> = T> + + X?_ are 
a w {t>J) = —NK j pj. 

T> + can be expressed in terms of D and N; so, V, T>, and iV contain all the information 
of K. It should therefore be possible to express the (ADM) gravitational Lagrangian 
(seeQ) 

= \ I N [R 2 - tr K 2 - R] e 

in terms of these variables. (Units are such that 4tvG = c = 1.) 

If we decompose the Lagrangian into "kinetic" and "potential" terms, then (^) 
immediately gives the potential term (n is again the dimension of E), 

L pot = -\ j^NRe= ^- 2 j N\V\ 2 ~\ 

The simplest expressions in Z>_ and T> yield the pieces needed to construct the 
kinetic term: 

<a pr (X> 2 )> = (N 2 k>p,k k p k ) = ±N 2 (K3k ) = l -N 2 tvK 2 

and 

(a pr (TO>_) 2 > = {{pNfrp,) 2 ) 

= ^+2) N2 ( ai Kja k k k + a 3 k k a k K 3 + a 3 K k a d k k ) 
= WW) N2 ( 2R2 + ^K k K ]k - a^k k k k a 3 ) 
= ^T)N 2 [2K 2 + {2-n)t,K 2 ]. 
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Combining these gives 



N 2 (K 2 -tiK 



so 



n(n+2) ) 2 + "( n ~ 4 ) 2? 2- P 2 



However, we want this in terms of T> not Z>_. If we replace P_ by T> = T>- + D + 
in the expression for Lk in then this can be expanded into a term quadratic in T>_, a 
term mixed in T> + and T>_ , and a term quadratic in T> + . The term quadratic in T>_ is 
simply the expression for L^ Q that we started with. In the mixed term, the principal 
symbol of the integrand contains precisely one factor of a, therefore its trace is and 
the term vanishes. For the term in T> + there is useful freedom to rearrange factors; 
since a pi (T> + ) is a scalar, we can commute T> + past anything else in this Wodzicki 
residue; so the term is 



n(n— 2) 



-n-2 



£>4 



Combining these results gives 



L 



N~ l \V\- n - A 



n(n+2) {^Y^jy^ _|_ n(n-i) jy2-p2 



I N- 1 V + \V\- n - 2 V + + ^ 



N\V 



2-n 



We still need to eliminate T> + from the second to last term of this expression. Inserting 
the definition (KJ) of T> + into that term gives 



n(n-2) 
16 



j N^p 2 , N]_\V\- n ~ 2 [V 2 , N]. 
j N~ l V 2 N\V\~ n - 2 [V 2 , N}_ 



n(n~2) 
16 



by expanding the first commutator. The last term vanishes because the principal 
symbol involved is an odd function of p. Expanding the remaining commutator gives 
directly 



16 



"( n ~ 2 ) J, N\V\ 2 ~ n - n(n ~ 2) i N-^NIVI^^NV 2 . 
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The final result is 



L 



N~ L V 



1 lTil-n-4 



n(n+2) ^~P) 2 + n ( n ~ 4 ) 2? 2 1? 2 



i / 6 , n(n-2) 
1 n-2 16 



jy|£>| 2 ~ n _ "(^ 2 ) NV 2 N- 1 V 2 N\V\~ i! ~ 2 



In the particular case of n = 3 this becomes 



iV- 1 !©]- 7 [f (VV) 2 - §D 2 X> 2 ] + f|^f - NV 2 N- l V 2 N\V\-\ 



In the case of n — 2, general relativity simplifies considerably. Unfortunately this 
formula for L cannot be applied in that case since n — 2 occurs in a denominator. 

Instead of calculating these Wodzicki residues geometrically, we can calculate them 
algebraically via the Dixmier trace. This allows us to apply this action to the general- 
ized case in which the spectral triple is no longer commutative. However, the first two 
terms here are written asymmetrically and so are not manifestly real in the noncom- 
mutative case. Fortunately, there is still some freedom to reorder the factors. One 
of the axioms of noncommutative geometry states that for any / 6 A, the operators 
[\V\, f]_ and [V, /]_ are bounded. This means that [/, |X>| _1 ]_ = \V\- 1 [\V\, f]_ {V^ 1 
is of order —2 — lower order than /|Z>| -1 . Since the first two integrands of (H) are in 
we can effectively commute iV _1 and |£>| _1 (or T>) within these expressions. 
Because of this, (|J) is real valued in general. 



8 Spatial Diffeomorphisms 

One striking difference between this form of Hamiltonian gravity and the ADM form is 
the absence of a shift vector. In the ADM formulation, the shift is needed because it is 
necessary to identify points of successive leaves of the foliation; the shift parameterises 
the freedom to choose that identification. In this algebraic formulation there is no 
need to refer to points; instead, spinor-functions are naturally identified on successive 
leaves using the Dirac equation. 

However, the shift plays another important role in the ADM formulation; variation 
of the action with respect to the shift gives the vector constraint - - the part of 
Einstein's equations concerning Goj. Since there is no shift in the algebraic action, 
it is not readily apparent that the vector constraint will be among the equations of 
motion coming from the algebraic action. It is, however, inevitable that the vector 
constraint will be present since in the constraint algebra it is generated from the 
scalar constraint given by variation of the lapse. It will be explicitly shown in this 
section how the vector constraint really does come from the variation of this action. 

The shift is a vector field, which is an infinitesimal diffeomorphism; so, in order to 
address this issue, it is necessary to first determine what a diffeomorphism is in the 
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algebraic formulation. The group of diffeomorphisms of a surface E is precisely the 
group of (outer) automorphisms of the algebra of functions C°°(£), but this group 
does not act naturally on the Hilbert space of spinors. Instead we need to use a 
group that does — namely, the group of automorphisms of the spinor bundle which 
preserve the density valued inner product. This is precisely the group of unitary 
operators which preserve the algebra of functions, the real structure, and (for n even) 
the grading. This group is essentially the semidirect product of diffeomorphisms with 
local unitary transformations of the spinor bundle. 

The generalization of a shift should live in the corresponding Lie algebra. In a 
small local region, the spinor bundle can be trivialized and spinor-functions look like 
half-densities valued in this bundle. A Lie algebra element looks like a Lie derivative 
plus an anti-Hermitian matrix-valued function, but the Lie derivative of a half-density 
with respect to £ is Ag := + • £) = |[£ fe , D k } + . The Lie algebra is 

flt = {He Op(ft) | E* = —E, [E,A t ]- C At, Ei = fCE, EJ = JE} 
= { At + X | £ G C°°(S t , T£ t ), X = -*t e C °°(J: t , End(SE t )), 

jx = xj, ix = xi} . 

In order to reintroduce degrees of freedom corresponding to a shift, drop the 
natural identification of Tit with Ti. and instead let the identification be variable. 
This is equivalent to composing the natural identification with a unitary operator 
U(t) from the above group; which is equivalent to replacing the fundamental variables 
in the action with U{t)V{t)U{t)- 1 and U{t)N{t)U{t)- 1 . Since the Lagrangian (§) is a 
Wodzicki residue of an expression in V, T>, and N, it is invariant under a simultaneous 
unitary transformation of these variables (treating T> as independent). The time- 
dependent unitary transformation of V and is thus equivalent to substituting V \— > 
V + [U^U.V]^ in the action. 

Now, start from U(t) = 1 (the natural identification) and infinitesimally vary U 
by H; this changes U~ l U by 5. This is thus formally equivalent to changing T> by 
[S, D]_ in the action Ldt. In order to extract the corresponding equations of 

motion, this should be expressed in terms of E not H; so, integrate by parts. This 
separates the change in the action into a time integral and the difference of a boundary 
term between the final and initial times. This boundary term is the change in the 
Lagrangian when T> is varied by [E,X>]_. 

In particular, if E = then the corresponding variation in T> is 

5(V) = [A^V]_ = ^[[e,D k ] + ,D J }^ 

= \a? [fc*,iy_,I>*] + + & [e,[D 1 ,Dk]-} + . 
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The last term simplifies, 

ia j [e, [Dj,D k ]J\ + = ia^ k [Dj,D k ] + i a j [[Dj, D k ]^,^ k ] _ 

= 

and so [A 6 £>]_ = -Jo* [[D j ,£ k ]„,D k ] + = -\ [aP^j,D k ] . Comparing this with 

the formula @ for T>, the effect on the Lagrangian is equivalent to varying K k - by 
N^ 1 ^^. This gives the change in the Lagrangian as 

which is exactly the vector constraint. Variation with H = X simply gives 0; this is 
in some sense an identically satisfied "Gauss' law" constraint. 

The action given by (§) doesn't actually know about the natural identification of 
Tit with 7i. Indeed, an infinitesimal unitary transformation of T> and N (coming from 
Q t ) is among the valid variations of these fields anyway. This means that the vector 
constraint will be among the equations of motion coming from the algebraic action, 
provided the variation is understood in the following sense. Normally an action 

L dt is required to be stationary with respect to variations with the arguments 
held fixed at the endpoints {to,ti} of integration; in this case the geometry at a 
given time is unchanged by a unitary transformation. If the variation is performed 
with the arguments at to and t% fixed only up to unitary equivalence, then there will 
be boundary-term equations of motion in addition to the Euler-Lagrange equations. 
These give precisely the vector constraint. This does also work in a more conventional 
setting; if the shift is simply set to in the ADM action, then the vector constraint can 
still be obtained by only fixing the initial and final geometries modulo diffeomorphism 
equivalence. 

9 The Noncommutative Generalization 

If the variables occuring in the algebraic action are restricted to come from a classical 
commutative space-time, then by construction this is completely equivalent to general 
relativity. However, the algebraic Lagrangian (H) can be applied equally well to a 
more general noncommutative geometry. In order to do this it is necessary to express 
the restrictions on the variables appearing in the action in a purely algebraic way, 
and to understand how the other structures of the spectral triple evolve with time. 
These other structures (At, J, and sometimes fi) enter into the action indirectly by 
restricting T> and N. 

The axioms of noncommutative geometry are actually not strong enough that in 
the ordinary manifold case the Dirac operator will be restricted to the form . We 
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actually have the freedom to add a (somewhat restricted) matrix-valued function to 
the Dirac operator without violating the axioms. This change is, however, detected 
by the Euclidean gravitational action; the matrix makes an ultralocal, positive defi- 
nite, quadratic contribution to the action. This means that the desired type of Dirac 
operator is distinguished within its equivalence class as minimising the Euclidean 
gravitational action. In the case at hand, the contribution to the gravitational La- 
grangian is the same ultralocal term multiplied by N. This means that the "correct" 
Dirac operators will again minimise the gravitational action, provided that the lapse 
is stricly positive. The equations of motion given by the algebraic action will contain 
the desired restriction on the Dirac operator, in addition to the Einstein equations. 

In order to understand how operators evolve in time, recall that this is the Heisen- 
berg representation; operators should evolve by the equation A = i[H,A]_. The 
Hamiltonian is determined by looking back to the Schrodinger representation. Sup- 
pose that a space-time spinor-function ift is time-constant with respect to some iden- 
tification of spinor bundles with shift. Then V n ip = Xip for some matrix- valued 
function X (independent of ip). The time derivative in the Heisenberg sense is given 
by the Hamiltonian, 

Hip = -if t ip = -iNifi) = -iN(V n - oPVj)tj) 
= [-iN(X - \K) + NV] i). 

This, together with the requirement that H must be self-adjoint, gives H = j[N, T>} + . 
It is trivial to add in a generalized shift S(t) G g t ', this gives the general Hamiltonian, 

H=±[N,V] + -iZ. 

The evolution of a specific f E At depends upon the choice of S, but since [S, a] _ G 
At, the evolution of At as a whole is independent of S. The evolution of At can be 
described axiomatically by the condition 

f(t) g \ [[N, V] + J]_ + A t V/ G C°° (R, AO • 

The S is irrelevant to the evolution of the other structures as well; for even di- 
mensions [S, fi]_ = implies that ft = | [[N, T>] + , if\_ = i[N,V] + yi. If N is strictly 
positive, then this formula can be inverted and yi(t) actually determines V. In the 
same way, for the n = 3 mod 4 case J = i[N, T>] + J, and J(t) determines T>. 

I have used the symbol C°° rather freely without really explaining what it means 
here. The most obvious definition of smoothness is the one associated with the norm 
topology: A time- dependent operator is C°° if all its derivatives are bounded operators. 
In the commutative case, the first derivative of a smooth section of Ar is indeed 
bounded, but the second derivative is generically unbounded. This definition fails for 
fl(t) where the first derivative is unbounded and for T> t which is itself unbounded. It 
is therefore necessary to use a much weaker definition of smoothness. 
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For any value of t we can define the dense subspace 7^°° C 7i as the common 
domain of all powers of T> t . In the commutative case Ti 00 is the space of smooth 
spinor-functions at E*, but since smooth initial data evolves smoothly under the Dirac 
equation this should be independent of t. By any reasonable definition, a smooth 
space-time operator acting on a smooth spinor-function should give a smooth spinor- 
function; this suggests the very weak definition for smoothness of a time-dependent 
operator: A G C°°(R,Op(H)) if and only if (i>\A\ip) G C°°(R) V*ij; 7 (p G H°°. If 
T> t satisfies this condition then 7i°° will indeed be independent of the t at which it 
is defined. This definition of smoothness is associated to a topology defined by the 
system of seminorms {^4 i— > |(t/>|A|t/>)| | ip G 7Y°°}. Since this is a subset of the 
seminorms which define the weak topology, this topology is even weaker than the 
weak topology. 

N is restricted by the assumption that C is time-independent. C = implies 
that C commutes with iH. Since C always anticommutes with i and T>, this means 
that iV must commute with C. Assuming that iV continues to be ^-even, this gives 
the condition N = CNC~ X = JNJ~ l . It is also reasonable to assume that N will 
be self-adjoint in general. In the commutative case the class of functions which 
commute with J is the same as the class of functions which are self-adjoint, namely 
real valued functions. Because of this there are many choices for the range of N 
which are equivalent in the commutative case. N(t) can be taken as a self-adjoint 
element that commutes with J in the algebra generated by At and JAtJ~ x , or in 
the sum A t + JA t J' 1 , or simply in At- The last, most restrictive case implies that 
N(t) G A t fl JA t J~ l C Z{At) (the center of At) since by the axioms JA t J" 1 is 
required to commute with At', this choice is actually supported by the special case of 
a Connes-Lott model. 

In a Connes-Lott type model of particle physics (see [|]]), the algebra is "slightly 
noncommutative" . It is a finite dimensional matrix algebra over the algebra of func- 
tions on a commutative manifold. The gauge group is a subgroup of the group of 
unitary elements in the algebra. Most choices of N would break this gauge symme- 
try, but if N is chosen from N(t) G A t to satisfy N = N* = JNJ~ l , then it will be 
proportional to the unit matrix — the space-time geometry will be simply the prod- 
uct of a commutative space-time with an "internal" geometry. On the other hand, 
with a broader choice of N, there might be interesting models possible in which is 
responsible for some symmetry breaking. 

In the commutative case the real structure does not allow the Dirac operator to be 
changed by a gauge potential. However, in the slightly noncommutative case a gauge 
potential can be added, and with this choice of iV it will not affect the gravitational 
action. This means that if the gravitational action is not supplemented with a gauge 
action then the dynamics of the gauge potential will be indeterminate. 

This example implies that in the general noncommutative case there will be per- 
turbations to which the action is insensitive, but this is not necessarily a problem and 
the precise significance can only be determined by looking at specific examples. 
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10 Time Diffeomorphisms 



In the ADM formalism there are many possible choices of foliation for any given 
space-time, and since they are all equivalent descriptions it is possible to transform 
from one foliation to another. 

An infinitesimal transformation of the foliation is determined by a real function 
77 G C°°(M.) or equivalently rj{t) G At- At each point of a time slice, 77 is the infinitesi- 
mal coordinate time that the slice is deformed by there; i]N is the infinitesimal proper 
time. The argument to determine the effect of this transformation on a space-time 
spinor-f unction i/j is exactly the same as the argument in section |2| for evaluating 
the time derivative, but with N replaced by rjN. This means that ifi transforms by 
Sip(t) = r)(t)ip(t). This determines the effect on any space-time operator A to be 
5 A = [77^, If / G C°°(1R, Ar) is a space-time function then it should transform 
by Sf = [rj4r, f]- — rjf + [77, f)-4;- In the commutative case the second term van- 
ishes identically and this is fine, but in the noncommutative case the second term 
will not vanish in general. A space-time function should transform to a space-time 
function in the new foliation; it cannot involve any time derivative operators. The 
only transformations of the foliation which are acceptable in general are those for 
which i](t) G Z(At)- This does, however, always include the reparameterisations of t 
which do not change the foliation. 

This accounts for the ugliness of the expression (^) for the gravitational action. 
In the ADM formalism there is a simple choice of variables such that N occurs only 
linearly in the action. As a result, the constraint equation that results from varying 
iV does not contain N itself, and it is possible to solve the equations of motion for 
the other variables with an essentially arbitrary choice of N. The arbitrariness of N 
corresponds directly to the arbitrariness of the foliation. In the general noncommu- 
tative case this is not so; N is not arbitrary so it must satisfy nontrivial constraints 
and it is inevitable that iV occurs nonlinearly in the action. 

If the most restrictive condition on iV is assumed, then the restriction on 77 should 
not be suprising; indeed, 77 should be restricted to the same form as N. If this 
restriction is chosen, then we should always be able to transform to iV = 1, at 
least for some finite time interval. This means that an arbitrary solution can still 
be obtained after setting N = 1 in the equations of motion (which simplifies them 
considerably). Unfortunately, as in the ADM formalism, an N = 1 description will 
generally only work for a finite time, since the time slices tend to develop cusps. 

11 Conclusions 

This model has rather mixed qualities. It inherits several good qualities from non- 
commutative geometry. Diffeomorphism invariance takes on a much simpler and more 
linear form than in ordinary differential geometry. This formalism eliminates the arbi- 
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trary identification of points which is necessary in other forms of Hamiltonian gravity, 
and with it the shift vector. 

On the other hand, the expression for the Lagrangian @ is ugly. In the passage 
from the classical case to the noncommutative generalization there are many poten- 
tially arbitrary choices; a commutator could in principle be randomly inserted almost 
anywhere. However, this may not be a problem with only T>, T>, and N appearing in 
the action. 

This model may be most useful as a step towards a toy model of regularized 
quantum gravity. If a field theory can be constructed on a noncommutative geometry 
with finite dimensional Hilbert space and algebra, then the theory cannot be divergent 
when it is quantized (see [0]). If the action can be written as the Dixmier trace of 
an operator in C^ 1,00 \ then this can be approximated (up to normalization) by the 
ordinary trace in the finite degree-of-freedom case; this is a direct result of equation 
(§). This is not the case for the gravitational action; the Dixmier trace, instead, picks 
out a logarithmically divergent part of a series which is faster than logarithmically 
divergent. 

The badly behaved part of the gravitational action is essentially just r Tr u) (\T>\ 2 ~ n N). 
This is calculated from the zeta function Cv,n(s). Suppose that the zeta function has 
only simple poles at s = n — 2 and s = n and has no other singularities to the right 
of s = n — 2; this is always the case for a commutative manifold. (The singularity 
set of the zeta function is called the "dimension spectrum" in f5j.) The residue of 
the pole at s = n is proportional to N, the space average of N. This means that 
Cv,(n-n)( s ) — (v,n(s) — N(d(s) is regular to the right of s = n — 2. This means 
that the trace of \V\ 2 ~ n (N — N) is only logarithmically divergent, at least in some 
averaged sense. So, Tr w |'D| 2 ~ n (which is simply the Euclidean gravitational action) is 
entirely responsible for the excessive divergence. If an appropriate approximation for 
the Euclidean gravitational action in the finite degree of freedom case can be found, 
then this probably could be easily extended to adapt (|8|). Although the Chamseddine- 
Connes spectral action H is done in that spirit, that formula cannot be used directly 
in this case since it approximates a higher derivative gravity, rather than pure general 
relativity. 

In summary, the geometry in this model is a generalization of space-time which 
consists of a smoothly time- dependent real spectral triple with a fixed Hilbert space, 
and a time-dependent "lapse operator" which comes from the algebra of functions. 
The real structure is time-independent for spatial dimensions n = 1 mod 4; for 
n = 3 mod 4 it is time-dependent; and for even dimensions it is time-dependent, 
but the product of the grading and the real structure is time-independent. The time 
dependences of the grading and real structures are determined by the Dirac operator 
and lapse. 

The gravitational action is given by a messy, second order (in j|) Lagrangian 
depending only on the Dirac operator, its time derivative, and the lapse. Although 
a first order Hamiltonian description exists in principle, the momentum conjugate to 
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the Dirac operator is valued in a difficult to parameterise space of linear functions of 
operators. The gravitational Hamiltonian probably cannot be written in closed form. 

In principle, in the cases when the grading or real structure are time-dependent, 
the Dirac operator is determined by this time derivative. This suggests that V could 
be eliminated in favor of J as a fundamental variable. However, this leads to a very 
nasty Lagrangian, quadratic in J, which cannot be written in closed form. That is 
probably not the appropriate way to go since J and fi do not appear directly in the 
original action. Instead, choose initial values for J and yi and evolve them using 
l 2 [N,V} + . 

In the commutative case it is indeed possible to reconstruct a description of the 
geometry in terms of the ADM variables. At each time t the point set topology 
is recovered from At, and the spatial metric can be reconstructed from T> t . If we 
choose a specific generalized shift S, then evolution by the corresponding Hamiltonian 
H = ^[N, T>] + — iE determines an identification between all the AtS. Identifying the 
functions is equivalent to identifying points, so this makes the geometry into a time- 
dependent Riemannian geometry on a fixed manifold. The lapse remains the lapse; 
the shift vector, acting as a derivation of functions, is a \— > [S, o]_. 

Appendix: Normalizations 

The proportionality constants of Tr w f\D\~ n oc J fe and (a_ n \D\ 2 ~~ n ) oc R can most 
easily be determined by looking at unit spheres. The sphere S n is an SO(n + 1) 
symmetric space; its spinor bundle is inherited from K n+1 and is therefore trivial. 
As a representation of SO(n + 1), the sections of the spinor bundle are the tensor 
product C[S n ] <g> S of the scalar functions C[S n ] with the (2 Li J -dimensional) spinor 
representation S; the scalar functions are the direct sum of all reduced symmetric 
powers of the fundamental representation. 

If we choose an origin o G S n , then there is a natural decomposition so(n + 1) = 
so(n) © t where the so(n) is the stabilizer of o, and t = T D S n is its orthogonal 
complement. Choose a basis {J a} C so(n + 1) of Hermitian generators that splits 
into {Jq} C so(n) and {Jk} C 6. Since this is a unit sphere, the Dirac vector at o is 
7fc = 2 (1 ® Jfc); this is fixed because j\j"fk] must give the spinor rotation matrices. A 
vector field £ G so(n+l) acts as a Lie derivative on spinor-functions, but the difference 
between a Lie derivative and a covariant derivative is something proportional to V£. 
The £ G t are precisely those £ G so(n + 1) such that V£ vanishes at o; therefore, 
the Lie derivative action of t at o is the same as the covariant derivative action. The 
covariant derivative at o is thus = — % Jk = —i (Jk <8> 1 + 1 <8> Jk) = —iJk <8> 1 — \ 7&- 
The Dirac operator at o is D = ij k Vk = 2Jk <E> J k + f , but since the action of J a 
on scalars vanishes at o, this is equivalent to D = 2 J a <8> J + f which is valid at all 
points. 

Using standard formulae for the quadratic Casimir (see JT2I), this gives the spec- 
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trum of the Dirac operator as Spec(-D) C Z + | with degeneracy 



(n-l) 



■2L5J 



|^|rt-l _ n(n-l)(n-2) |^|n-3 _|_ 



Note that the lowest eigenvalue of D 2 = \R - V 2 is > = as it 

should be. The zeta function is (the analytic continuation of) Cd(s) = Tr |-D| _S = 
J2\ [d(ty + ^( — ^)] — 2~2\ 2g?(A)A~ s , where the sum is over positive integers (half 
odd integers) for even (odd) dimension. If the degeneracy is written d(X) = ao|A| n_1 -|- 
ai|A| n ~ 3 + . . . , then the zeta function can be written 

Cd(s) = J2 2 M™- 1 -* + a!A"- 3 - s + ...]= 2a kC(s -n + 2k + l) 

A k=0 

where, for even dimension, £ is the Riemann zeta function. For odd dimension, replace 
((s) with (2 s — l)C(s); both of these functions are meromorphic, with a single simple 
pole at s — 1 of residue 1. This shows that Cd( s ) is a meromorphic function with 
simple poles at s = n, n — 2, . . . > 0; the residue at s = n — 2k is 2a^. 

The residue at n shows that TrJ_D|~" = ^ Res s= „ (d(s) = ^a = ^LfJ. Com- 
paring this with 

D\- n = / e 



,s,. (V)' 
shows that (in general) 



A - lZ_j K2> Tr A 

2 L 2 J 



Comparing the residues at n — 2 and n shows 

i|£)|2-n = gl ipj-n = _ n(n-l)(n-2) /* e = f Re , 

J a J JS n JS" 

This means that (in general) 
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